Effect of Particle-Hole Asymmetry on the Mott-Hubbard Metal-Insulator Transition 
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The Mott-Hubbard metal-insulator transition is one of the most important problems in correlated 
electron systems. In the past decade, much progress has been made on examining a particle-hole 
symmetric form of the transition in the Hubbard model with dynamical mean field theory where it 
was found that the electronic self energy develops a pole at the transition. We examine the particle- 
hole asymmetric metal-insulator transition in the Falicov-Kimball model, and find that a number of 
features change when the noninteracting density of states has a finite bandwidth. 



PACS numbers: 71.30.+h, 71.27.+a, 71.10.Fd 

The Mott-Hubbard metal-insulator transition [1,2] is 
a classic example of the physics of strongly correlated 
electrons. The physical mechanism of the transition 
arises from a local Coulomb repulsion U that forbids 
double occupancy of the electrons, creating an insulator 
when there is one particle per site (on average). Exper- 
imentally it is found in a variety of materials, includ- 
ing many transition metal compounds (MnO, NiO, NiS, 
YBa 2 Cu306, etc.) for which band structure calculations 
severely underestimate the band gap or yield a metallic 
solution [3]. 

In the Hubbard model [2] , a number of approximations 
have been employed that assume either the metallic state 
and head toward the insulator [2,4,5], or vice versa [2,6]. 
Often it is difficult to develop an approximate theory 
that is able to describe both the weakly correlated Fermi- 
liquid phase and the strongly correlated insulator, hence 
most approximate methods yield only limited informa- 
tion about the transition. Much progress has been made 
with dynamical mean field theory (DMFT) [7], where 
fundamental questions such as, is the transition contin- 
uous or discontinuous, does the Fermi-liquid metal sur- 
vive up to the transition or do non-Fermi-liquid metallic 
phases intervene, do metastable phases exist, and so on 
have been analyzed in great detail. However, the nu- 
merics are quite complicated and delicate, because the 
different phases are separated by very small energies. 
Hence there has been much controversy about the an- 
swers to these questions and about the details of the 
Mott-Hubbard transition [8-15]. 

Some of these problems arise from the limit of infinite 
dimensions, where the noninteracting density of states 
(DOS) on a hypercubic lattice is a Gaussian, which has 
an infinite bandwidth. This means that the DOS can 
only vanish (and thereby yield an insulator at T = 0) 
when the self energy diverges. Since this occurs only at 
the single point of an isolated pole, the DOS in the insu- 
lating phase is really a pseudogap, and is nonzero (albeit 
exponentially small) in a region around the pseudogap. 
Hence, the MIT on a hypercubic lattice always occurs 
when the self energy develops a pole (that lies infinites- 
imally below the real axis). It turns out that the same 



scenario occurs on the Bethe lattice at half filling, even 
though the DOS has a finite bandwidth the MIT oc- 
curs when the self energy develops a pole at the chemical 
potential (although here there is now a well-defined gap, 
where the DOS vanishes over a finite range of frequen- 
cies). 

The question we wish to address is what happens 
when the MIT occurs for a system that does not pos- 
sess particle- hole symmetry [16]. Following the argument 
given above, the pole formation must be the underlying 
cause of the MIT on the hypercubic lattice, but the situa- 
tion on the Bethe lattice is unclear. The interest in exam- 
ining particle-hole asymmetric cases lies in the fact that 
most real materials do not have particle-hole symmetry, 
so understanding consequences of breaking particle-hole 
symmetry is important for understanding experimental 
systems. There are two ways to break particle-hole sym- 
metry: (i) one can modify the lattice, so it is not bi- 
partite, then there is no particle-hole symmetry at half 
filling, where the MIT occurs, or (ii) one can modify the 
model so that a MIT occurs at different fillings, away 
from the particle- hole symmetric limit. We choose to 
examine the latter here. 

We consider the MIT in the spinless Falicov-Kimball 
(or simplified Hubbard) model [17], which is believed to 
describe correlated-electron behavior, and the MIT in 
particular, in materials that can be fit into a binary alloy 
picture. The canonical system that fits this picture is 
Ta^N [18,19] which has its MIT occur at a particle- hole 
asymmetric value of x = 0.6. The Falicov-Kimball model 
has been previously applied to the MIT [20,21], its ad- 
vantages over the Hubbard model are that, it has a MIT 
for a wide range of fillings and the numerics are under 
much better control. The Falicov-Kimball Hamiltonian 
has the following form 



(1) 



where cj (q) denotes the creation (annihilation) opera- 
tor for a spinless electron on site i, and the summation is 
restricted to nearest neighbors. The classical variable Wi 
equals one or zero, corresponding to the presence of an 
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A or B ion at the given lattice site, and U is the diagonal 
site-energy difference between the two ionic configura- 
tions. The hopping integral t is appropriately scaled in 
order to be finite in the limit of large dimensions (d) or 
coordination number (Z) [7], and to result in the same 

effective bandwidths W = 



J e 2 p(e)de (with p(e) the 

noninteracting DOS) for both lattices [14]. Therefore, 
t = t*/s/~Z on the Bethe lattice and t = t*/V2d on 
the hypercubic lattice. The conduction electrons inter- 
act with the localized particles (which have an average 
filling of w\ = (wi)) with an interaction strength U; this 
forms the canonical binary alloy picture. Such a picture 
is particularly useful because the model exhibits a MIT 
when the particle-hole symmetry is broken (w\ =/= 0.5), 
as long as the total particle density equals one. Hence, 
we constrain our calculations to fix the total number of 
particles, i.e. the number of conduction electrons satis- 
fies p e = 1 — w\. If we choose to measure all energies 
in the units of t* = 1, then the noninteracting DOS on 
the Bethe and hypercubic lattices are a semicircle and 
a Gaussian, respectively, PBethe(t) = V4 — e 2 /2ir and 
p H c(e) = exp(-e 2 /2)/ v / 2^. 

In the limit of infinite dimensions the Falicov-Kimball 
model can be solved exactly using DMFT [22,23]. Be- 
cause the self-energy S(w) has no momentum dependence 
one can employ an iterative scheme using the following 
relationships between the self-energy S(w), the retarded 
Green's function G(u), and the effective medium Gq(w) 
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closed in the iterative loop until the self-energy is con- 
verged to a desired accuracy [24]. The algorithm yields 
reliable results for relative accuracies up to 10~ 13 . The 
calculations are particularly easy in the case of the Bethe 
lattice since the integral in Eq. (2) can be determined an- 
alytically (for the hypercubic lattice it is a complex error 
function). Once the algorithm is converged, the interact- 
ing DOS is defined to be pmt{w) = —lm[G(u))]/n. The 
essential difference in the physics on these two lattices 
arises from the behavior of the interacting DOS. While 
the Bethe lattice exhibits a well defined gap where the 
interacting DOS is exactly zero for a finite range of fre- 
quencies when U > U cg (the critical value of the interac- 
tion strength for opening a gap in the interacting DOS), 
on the hypercubic lattice the DOS splits into two sub- 
bands divided by a "pseudo-gap" , where the interacting 
DOS is exponentially small and exactly zero only at one 
point. The interacting DOS away from half filling has 



appeared elsewhere [20,25]. We should note that in the 
Falicov-Kimball model there is no quasiparticle peak at 
the chemical potential, which normally develops for U 
less than U cg in the Hubbard model in infinite dimen- 
sions. This is because the Falicov-Kimball model lacks 
quasiparticles in the strict definition of the term, since 
the lifetime of Fcrmionic excitations never becomes infi- 
nite as T —* 0. 

The metal-insulator transition at half-filling is closely 
related to the development of a pole in the self-energy, 
where the DOS is suppressed to zero (U cg ) at the same 
interaction strength for which the pole forms (U cp ); i.e., 
U C g = U cp . Therefore one could suggest to use the residue 
of the pole as an order parameter for the MIT. In fact, 
a plot of the residue versus (U — U C p)/U C p is universal, 
for all fillings on both the hypercubic and Bethe lattices, 
indicating that a scaling theory holds for the residue of 
the pole (see the inset to Fig. 1). Away from half-filling, 
these two processes (pole formation and the MIT) are 
decoupled on the Bethe lattice, with the pole formation 
occurring after the MIT (in particular, the real part of 
the Green's function does not cross the horizontal axis 
[within the band gap] until the pole forms). 

We begin with the cubic equation which is satisfied by 
the Green's function on the Bethe lattice [20] 

G 3 -2a;G 2 + (l + a; 2 -^)G-(a ; + a) = (5) 

where x = u + p — U/2 and a = U{w\ — \). The pole in 
the self-energy develops when G = 0, requiring x + a = 0, 
or ujpoie = U(l — w\) — p. Next, one notes that G = is 
a physical root [26] when 1 + x 2 — [7 2 /4 < 0, which leads 
to the simple formula for U cp , the critical value of U at 
which the self-energy develops a pole 
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cp 



(6) 



Away from half filling, the critical values of U for the 
MIT and the pole formation are not necessarily the same 
on the Bethe lattice, and the former is found from the 
following. Let us use standard notation for the cubic 
equation's coefficients (see [27]) q = (3 — x — 3C/ 2 /4)/9 
and r = x(l + U 2 /2)/6 + a/2 - x 3 /27. The condition 



for the location of the band edges (q 3 



0) is a 



fourth order equation in u, and must have exactly three 
distinct roots at U cg the critical interaction strength for 
the MIT (four distinct roots for U > U cg ). As shown 
in Ref. [20] this condition leads to an equation for the 
critical interaction: 
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a 2 4(U 2 - 1) 



U 2 



27U 2 



(7) 



Solving Eq. (7) for U yields the critical interaction 
strength for the MIT (gap opening on a Bethe lattice) 
U, 



eg ■ 



2 



U cg = yjl + 3w\ /3 (l - Wl )V3[(l - Wl )i/3 + w \'\ (8) 

Both expressions (6) and (8) yield 2 at half filling (wi = 
0.5). The behavior of U cp and U cg for arbitrary fillings is 
shown in Fig. 1. 
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FIG. 1. Phase diagram for the MIT in the Falicov-Kimball 
model on the Bethe and hypercubic lattices. The solid line 
denotes U cg for the gap formation on the Bethe lattice, the 
dashed line denotes U cp for the pole formation on the Bethe 
lattice and U cg = U cp on the hypercubic lattice. The inset 
shows the universal curve of the residue of the pole as a func- 
tion of the interaction strength R(U) for both lattices and all 
fillings. 

As the system moves away from half-filling the particle- 
hole asymmetry allows for the formation of a third 
"phase", in which the interacting DOS has a gap but 
there is no pole in the self energy (region between the 
solid and dashed lines in Fig. 1 on the Bethe lattice). 
This dramatic difference in the MIT as we move away 
from half filling (where the transition occurs at smaller 
U values than that at half filling on the Bethe lattice, but 
larger U values on the hypercubic lattice) is likely tied to 
the fact that there is an infinite number of neighbors and 
loops on a hypercubic lattice, which allows one to avoid 
sites blocked by the presence of a static particle, but on 
the Bethe lattice, whole subtrees become disconnected 
because there are no loops, so the blockage is more effec- 
tive. On the hypercubic lattice, the pole formation and 
MIT always occur at the same value of U (U cg — U cp ) 1 
and it is possible to obtain an exact expression for the 
pole formation, which turns out to coincide with Eq. (6). 
Numerical calculations show that the interacting DOS for 
the hypercubic lattice has an exponentially small region 
over a similar range in frequency as the gap region on the 
Bethe lattice (for the same value of U). 

Fig. 2 shows the relative interaction strength as a func- 
tion of the relative location of the pole within the gap, 
on the Bethe lattice. The pole is located in the middle of 
the gap at half filling for any value of U > U cg - However, 



as the particle-hole symmetry is broken, the pole first 
appears at the lower or upper band edge (for U = U cp ), 
depending on whether w\ is larger or smaller than 0.5, 
and as U increases the pole drifts closer to the center of 
the bandgap. Note that there is no smooth transition 
between the half-filled case and the particle-hole asym- 
metric case. 
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FIG. 2. The relative interaction strength for the opening 
of a gap versus the relative location of the pole (E po i e ) in the 
gap (of width E gap ) on the Bethe lattice; the origin on the 
horizontal axis lies at the center of the gap. The lines from 
left to right correspond to different fillings toi ranging from 
0.9 to 0.1, in steps of 0.1. 

Numerical calculations show the evolution of the real 
part of self-energy on the Bethe lattice compared to that 
on the hypercubic lattice in Fig. 3, when approaching the 
MIT from the metallic side. Aside from obvious similari- 
ties between the two lattices, this graph demonstrates the 
existence of the third phase in the middle panel, where 
U = 2. The half-filled curves (wi = 0.5) always show 
poles in the insulator, and therefore, the DMFT scenario 
discussed in Refs. [5,14,11] holds. However, the curves 
corresponding to wi = 0.25 exhibit large (negative) but 
finite values of Re[E(o>)] (the pole has not yet developed) 
but the system is an insulator with a well developed gap 
at this value of U. Hence, we conclude that the develop- 
ment of the pole and the MIT are decoupled away from 
half filling and even if it might be tempting to use the 
residue of the pole as an "order parameter" for the MIT, 
it fails to describe the situation off of half filling on lat- 
tices with a finite bandwidth. So what significance can 
be made of the pole formation? In order to investigate 
this, we calculate the dc conductivity in the relaxation 
time formalism (see [28] for details). These calculations 
(not shown here) indicate that the MIT is always a con- 
tinuous transition at T = 0, with the conductivity being 
suppressed continuously to zero as U — > U cg . Also, at 
finite temperature we see no evidence for the influence of 
the pole on the transport. The conductivity curves are 
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smooth functions of T with no unusual features occurring 
when U > U cv . 



Bethe hypercubic 




FIG. 3. Evolution of the real part of self-energy for different 
values of U, for both the Bethe (left panel) and hypercubic 
(right panel) lattices, at half filling wi = 0.5 (solid line) and 
at wi — 0.25 (dashed line). The origin of the frequency axis 
is set to U(l — wi) — jU. 

We need to note that in the correlated insulator, in 
the vicinity of u> = £7(1 — wi) — fj,, one can obtain ana- 
lytic expressions for the Green's function G{ui) and the 
self-energy on the hypercubic lattice [25]. As ex- 

pected, one obtains imaginary parts of G(u>) and S(w) 
exponentially approaching zero as uj — > U(l — tui) — fi, 
while they are exactly zero on the Bethe lattice. These, 
at first sight, small differences, in fact turn out to cause 
dramatic changes in the low temperature transport prop- 
erties [25,29] when going from the Bethe lattice (where 
the relaxation time vanishes in the gap) to the hypercubic 
lattice (where the relaxation time has a power-law depen- 
dence not an exponentially decaying dependence around 
the pseudogap), and therefore must be dealt with care- 
fully The most significant departure is in thermal trans- 
port properties, where the thermopower diverges on the 
Bethe lattice, but vanishes on the hypercubic lattice in 
the insulating phase as T — > (when particle-hole sym- 
metry is broken). 

In conclusion, we have analyzed the effect of particle- 
hole asymmetry on the Mott transition in the infinite 
dimensional Falicov-Kimball model (on both the Bethe 
and the hypercubic lattices). Hitherto, it was believed 
that the scenario for the MIT on both lattices was the 
same, as indeed is the case at half-filling. We find that 
this is not true when the particle-hole symmetry is re- 
moved, as is often the case in real materials. We show 
that in the absence of particle-hole symmetry the pole 
formation and the MIT are two unrelated processes on 
the Bethe lattice. So even though the residue of the pole 
satisfies many of the properties expected of an order pa- 
rameter, it cannot be employed to describe the MIT in 



all cases (although one could use it on the hypercubic 
lattice). Furthermore, there seems to be little difference 
between the properties of a correlated insulator with or 
without a pole in the self energy. We conjecture that 
all of the conclusions about the character of the MIT on 
the Bethe lattice will hold for other systems with a finite 
bandwidth, and hence much of these results will play a 
role in realistic models of the MIT for real materials in 
finite dimensions. 
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